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We apply the generating function technique developed by Nazarov to the computation of the 
density of transmission eigenvalues for a two-dimensional free massless Dirac fermion, which, e.g., 
underlies theoretical descriptions of graphene. By modeling ideal leads attached to the sample 
as a conformal invariant boundary condition, we relate the generating function for the density of 
transmission eigenvalues to the twisted chiral partition functions of fermionic (c = 1) and bosonic 
(c = — 1) conformal field theories. We also discuss the scaling behavior of the ac Kubo conductivity 
and compare its different dc limits with results obtained from the Landauer conductance. Finally, 
we show that the disorder averaged Einstein conductivity is an analytic function of the disorder 
strength, with vanishing first-order correction, for a tight-binding model on the honeycomb lattice 
with weak real- valued and nearest-neighbor random hopping. 



I. INTRODUCTION 

The recent manufacture of a single atomic layer of 
graphite (graphene) has renewed interest in the trans- 
port properties of Dirac fermions propagating in two- 
dimensional spacei 2 -' 3 -'^ Recent theoretical work in- 
cludes, among many others, the computation of the Lan- 
dauer conductance for a single massless Dirac fermion 
by Katsnelson in Ref. [H, as well as the computation of 
the Landauer conductance and of the Fano factor in 
Ref. H confirming the result of Ref. d and predicting 
sub-Poissonian shot noise. 

The conductivity has long been known to be related 
to a twist of boundary conditions^ This idea has been 
further developed by Nazarov who proposed a generat- 
ing function for the density of transmission eigenvalues in 
quasi-one-dimensional disordered conductors 'li . With 
this formalism, Lamacraft, Simons, and Zirnbauer repro- 
duced in R ef. ITU ( see also Ref. [13) nonperturbative results 
of Refs. [MillB for the mean conductance and the den- 
sity of transmission eigenvalues of quasi-one-dimensional 
disordered quantum wires for three symmetry classes of 
Anderson localization. 

The purpose of this paper is to establish a connec- 
tion between (i) the density of transmission eigenvalues of 
the noninteracting Dirac Hamiltonian describing the free 
(ballistic) propagation of a relativistic massless electron 
in two-dimensional space and (ii) twisted chiral partition 
functions of a combination (tensor product) of two con- 
formal field theories (CFTs) with central charges c = 1 
and c = — 1. In this way we provide a complementary 
method for calculating the Landauer conductance of a 
single massless Dirac fermion which agrees with the di- 
rect calculations of Refs. iand@, while it might give us a 
powerful tool to account for the nonperturbative effects 
for certain types of disorder. This connection comes from 
the observation that ideal leads attached to the sample, 



which are necessary to define the Landauer conductance, 
can be replaced by a set of conformally invariant bound- 
ary conditions. That this is possible is a consequence 
of the "diffusive" nature of the ballistic Dirac transport, 
as epitomized by its sub-Poissonian shot noise, which in 
turn suggests an insensitivity to the modeling of ideal 
leads. 

The transport properties of noninteracting Dirac 
fermions in two-dimensional space under ideal condi- 
tions, i.e., without the breaking of translation invari- 
ance by disorder, have also been discussed in terms of 
the Kubo formula. However, extracting the dc value 
of the conductivity from the Kubo formula, which de- 
pends on frequency to, temperature 1//3, and smearing 
r\ (imaginary part of the self-energy), is rather subtle, 
as it is known that the conductivity in the dc limit 
(tj, 1/(3, 77) = (0,0,0) is sensitive to how one approaches 
the dc limit, 16 ' 17 ! 18 . 19 ! 20 . 21 ! 22 ! 23 ! 24 There have been at 
least two known limiting procedures: The Einstein con- 
ductivity, defined by first taking the zero temperature 
and then the zero frequency limits while keeping r\ finite, 
is given by 1/tt in units of e 2 / h .ItlLl&ia ^2^4 n 
the other hand, by switching the order of zero smearing 
and zero frequency limits, one obtains, instead, n/8 for 
the dc conductivity (again in units of e 2 //i)j 16 i 17 

To clarify the origin of this sensitivity to how the dc 
limit is taken, we compute the ac Kubo conductivity 
without taking any biased limit in u>, 1/(3, or 77. Due 
to the scale invariance of noninteracting Dirac fermions, 
we show the Kubo conductivity to be a scaling func- 
tion of two scaling variables. We discuss several limit- 
ing procedures, including the above two, and clarify the 
relationship between different dc values. In particular, 
we demonstrate that, if we take the zero smearing limit 
prior to the other two limits, we can obtain any value 
between and ir/8 for the dc Kubo conductivity. The 
Einstein conductivity agrees with the conductivity deter- 
mined from the Landauer conductance. We also compute 



2 



several asymptotic behaviors of the Kubo conductivity 
that have not been obtained before. These considera- 
tions may be of relevance to experiments on graphcne if 
different limiting procedures are accessed. 

The perturbative effects of disorder in the form of 
weak real- valued random (white-noise) hopping between 
nearest-neighbor sites of the honeycomb lattice at the 
band center are discussed. We show that, as a conse- 
quence of the fixed point theory discussed in Ref. Hf| the 
Einstein conductivity is an analytic function of the disor- 
der strength. We also show that the first-order correction 
to the Einstein conductivity vanishes, in agreement with 
a calculation performed by Ostrovsky et al. in Ref. [26|. 
This result is of relevance to the two-dimensional chiral- 
orthogonal universality class . 27 ] 28 i 29 i 30 Potential relation- 
ships of the two-dimensional chiral symmetry class with 
certain types of disorder in graphcne were recently dis- 
cussed in Ref.[26|. We refer the reader to Refs.[26|,[3l|, and 
[32I for a discussion of white-noise disorder in graphe ne in 
terms of symmetry classes and to Refs. [3lll32H33U34ll35l 
for the possibility that smooth disorder could induce 
crossovers between different symmetry classes. 



II. 



MODEL 



Our starting point is the single-species (or one-flavor) 
Dirac Hamiltonian 



JC := -ihvpa^d^, 



H := J d 2 r&'H$ , (2.1) 



where is a two-component fermionic creation 

(annihilation) operator and v F the Fermi velocity. Wc 
choose ct„_„ to be the first two of the three Pauli ma- 
trices a x , <7 y , and a z in the standard representation. (We 
use the summation convention over repeated indices.) 
Hamiltonian (|2.ip describes the free relativistic propa- 
gation of a spinlcss fcrmion in two-dimensional space 
parametrized by the coordinates r = (x,y). As such, 
it possesses the chiral symmetry 



(2.2) 



The single-particle retarded and advanced Green's 
functions are defined by 



5^ /A (e) := {e±ihr,-K)-\ 



(2.3) 



Consequently, at the band center e = 0, they are related 
to each other by the chiral transformation as 



a z S*(e = 0)a z 



0). 



(2.4) 



Below, matrix elements between eigenstates of the posi- 
tion operator of the single-particle retarded Green's func- 
tion evaluated at e are denoted by 9^ (r, r';s). 

In the presence of an electromagnetic vector potential 
Ap(r), one modifies Hamiltonian (|2.ip through the min- 
imal coupling d — > d — i(e/Hc)A (with e < 0). The 



conserved charge current then follows from taking the 
functional derivative with respect to A (r), 



5'/i( r ) = ( r )> U :z 



' J F <V 



(2.5) 



There is no diamagnetic contribution due to the linear 
dispersion. 



III. KUBO AND EINSTEIN CONDUCTIVITIES 

A. Linear response 

We start from the bilocal conductivity tensor at a finite 
temperature 1/(3 defined by the linear response relation 
in the frequency-w domain 

f» d (r,cj,(3)= J d 2 r' a^{r,r'- 1 u J ,5^)E u {r\u J ) (3.1a) 

between the u> component E v {r' ,uS) of an electric field 
that has been switched on adiabatically at t = — 00 and 
the induced local current j'!f d (r, w, 0), where^ 



(3.1b) 



Here, 



(3.1c) 

is the response function, j M (r, t) is the current opera- 
tor in the Heisenberg picture, and (• • • )g is the expecta- 
tion value taken with respect to the equilibrium density 
matrix at temperature 1//3. The small positive number 
5 > implements the adiabatic switch-on of the electric 
field. The conductivity tensor in a sample of linear size 
L is defined by integrating over the spatial coordinates 
of the bilocal conductivity tensor, 



1? 



d 2 r'a„ v (r,r';cj,5,p). (3.2) 



We impose periodic boundary conditions and choose 
to represent the Dirac Hamiltonian (|2.ip by 



p a— ± 



(3.3a) 



where the fermionic creation operators a£ CT with a = 
± create from the Fock vacuum |0) the single-particle 
eigenstates with momentum 



\m) = \p,a) = &IJ 



V2\P\ 



\Px - W v 



\P\ 



(3.3b) 
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of "K with the single-particle energy eigenvalues e 



functions are given bys£ 



av F \p\ 



av. 



Pi 



(3.3c) 



The conductivity tensor can then be expressed solely in 
terms of single-particle plane waves 

a fiv (uj,S,f3,L)=^ Jde J^J dV 

x ^{M^( r )\ n )( n \h( r ')\ m ) (3.4) 

//3( £ )-//3( £ n) 



x S(e-e m ) 



e — e n + H(uj + iS) ' 



where /^(e) := + 1) is the Fermi-Dirac function at 

temperature 1//3 and at zero chemical potential. 

As usual, the infinite-volume limit L — ► oo has to be 
taken before the S — > limit in Eq. (|3.4p . (Recall that 
(5 controls the adiabatic switching of the external field.) 
In the following, it is understood that we always take 
these limits prior to any other limits. We thus drop the 
explicit S and L dependence of the conductivity tensor 
henceforth. The dc conductivity can then be computed 
by taking the subsequent limit, uj — > 0. The temperature 
1/(3 can be fixed to some arbitrary value. 

The real part of Eq. (|3.4p can be further rewritten in 
terms of single-particle Green's functions. This can be 
done by first replacing the two 5 functions in the real 
part of Eq. (|3.4[) . which appear after taking the 6^0 
limit, by two Lorentzians with the same width hq (see, 
for example, Ref. H3). Then, each Lorentzian can be 
rewritten as the difference of the retarded and advanced 
Green's functions, 9 R {e) — Sjf(ff)- By also noting that 
the transverse components (/j, ^ v) of the conductivity 
tensor (|3.4p vanish by the spatial symmetries of the ma- 
trix elements in J d r J d 2 r' {m\j^(r)\n){n\j v (r')\rn), we 
obtain 



Reo- Mi/ (cj,77,/3) =5^ — I de 



fp(s + hui) - f p {e) 



where we have introduced 



^ v (r,r';E,cj,ri) :=tr [5 A ~ R (r', r; e)j M 

x g ? f->,r';e + Mj, 



(3.5) 



(3.6) 



and j is defined in Eqs. (|2.5p . Here the trace is taken 
over spinor indices and 

9^ R (r, r'; e) := S>, r'; e) - S R (r, r';e). (3.7) 

Using translational invariance, the single-particle Green's 



3 R/A (k;s) 



Hr 2 - ri ) S R/A ik . £ ^ 



(3-8) 



e ± ihr] — hv F k ■ a 



In order to define the conductivity in the clean system 
we should take the r\ — > limit before the uj — > limit. 
On the other hand, rj can be interpreted physically as a 
finite inverse life time (imaginary part of the self energy) 
induced by disorder. Thus, it is meaningful to discuss Eq. 
(|3.5[) in the presence of finite rj. Below, we first discuss 
the 77 — ^ limit. We will then discuss the case of finite 77. 



B. ? 7 -> limit 

We define the ac Kubo conductivity tensor at any finite 
frequency uj > and temperature by 



^(w,/?) ■= lim Re a (uj,t],(3). 



(3.9) 



With the help of Eq. (3J 

n(ev F ) 2 



^(m|CT /i |n)(n|<7 t 



X [//j(0 -//3( £ J] S i £ m ~£n 



(3.10) 



When uj > and (3 < 00 , the sum over the basis (13.31)1) 
in the real part of Eq. (|3.10|) can be performed once 
the matrix elements of the currents have been evaluated, 
yielding 



(3huj 



(3.11) 



Observe that Eq. (|3 . is independent of the Fermi ve- 
locity v F £2- Finally, the ac Kubo conductivity (|3.1ip de- 
pends solely on the combination 



(3hw G R. 



(3.12) 



The limiting value of Eq. p. lip when uj — * and (3 — > 
00 can be any number between and e 2 7r/(8/i) provided 
the scaling variable (|3 . 1 2[) is held fixed. For example, if 
the limit (3 — > 00 is taken before the limit uj — > 0, then 

hm lim o*(u,l3)=6 lu ,^l. (3.13) 

This limitin g p rocedure reproduces the results from 
Refs. HH and Il7l. On the other hand, if the limit uj — > 
is taken before the limit (3 — > 0, then 



lim lim^>,/J)=0. 



/3^oo oj->0 



(3.14) 
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Clearly, liraa^ of„(u>, /3) = S fll/ {ir/8)(e 2 /h) for any fi- 
nite frequency u>, while lim^^o cr Jf l /( w , /3) = for any 
finite temperature The singularity at lu = 1/(3 — 
is a manifestation of the linear dispersion of the massless 
Dirac spectrum leading to a dependence on the scaled 
variables ui and /3. 

C. Case of finite r\ > 



Whenever i] > 0, Eq. (|3.5[) is an analytic function of A 
and Y at (A, Y) = (0, 0). For X, Y < 1, the power series 
expansion of Eq. (|3.5[) in terms of X and F is given by 



Recr xx (X < 1,Y < 1) = 




(3.21) 

up to terms of order X 4 , Y 4 , or X 2 Y 2 . 



For 77 finite, it is shown in Appendix [X] that the real 
part of the longitudinal conductivity, Eq. (|3 . 5|) , is a scal- 
ing function of two variables, i.e., 



Rea xx (u J , V ,p) = Rea xx (X,Y), (3.15) 



where 



X 



Y 



1 

/3/V 



(3.16) 



1. dc response uj — 

If we take the dc limit uj — > while keeping 77 > 
finite, Eq. (|3 . 5[) can be expressed as 



Re^(X = 0,r)=^A fde^f^ 



x ^(rj/je.w = 0,77), 



where is defined in Eqs 
help of [see Eqs. (ED} and (j2~5]) ] 



(3.17) 
and (|3~7)) . With the 



j' M = ie [5£, /ft, 



one can show that^ 



Rea Mi/ (X = 0,Y)=5^- V 2 J de^ 

xtr [g^(0,r;e)g^(r,0;e)] 



j2 2 

a rr 



(3.18) 



(3.19) 



with fj, = x,y. Equation (|3.19|) is usually referred to as 
the Einstein conductivity since it is related to the diffu- 
sion constant via the Einstein relation (see, for example, 
Ref. HO). 

A closed-form expression for Eq. (|3.17| can be obtained 
at zero temperature, 

Be<T ltu (X = 0,Y = 0) = 5^~. (3.20) 

The same value was derived in Ref. [lfl Related predic- 
tions were also made, among others, in Refs. HH, H3, Hl|, 
S3, H, [M and [13^ Equation (|3~2^]) also agrees with 
the conductivity determined from the Landauer formula. 
(This was first observed in Ref. [f| We will reproduce 
this fact in Sec. IIVI from Nazarov's generating function 
technique.) 



2. Arbitrary u> and /3 at finite r\ 

For generic values of (X — u>/t],Y = l/(3hr]), we are 
unable to evaluate the real part of the longitudinal con- 
ductivity from Eq. (|3 - 5[) in closed form. The numerical 
integration of Eq. (|3.5[) when /j = v = x or, cquivalcntly, 
of Eq. (|A2a|) is presented in Fig. [Qi^ 

First, we fix Y = l/(pTvq) and discuss the X = uj/?] 
dependence of the real part of the longitudinal conduc- 
tivity from Eq. (|3.5p . Wc distinguish two limits. When 
X — > 00 as happens when the frequency uj is much larger 
than the inverse life-time 77, the real part of the lon- 
gitudinal conductivity measured in units of e 2 /h con- 
verges to the limiting value ir/8. In the opposite limit 
of A — > 0, the limiting value is an increasing function 
of Y and is given by Eq. (|3 . 2 1 [) with X = 0, when 
the temperature 1//3 is much smaller than the energy 
smearing fey {Y <C 1). These two limiting behaviors are 
smoothly connected as is illustrated in Fig.QJa). For ex- 
ample, (h/e 2 ) Rea xx (X, Y = 0) increases monotonically 
between 1/77 and ir/8 as a function of A^ The approach 
to the limiting value ir/8 for large A when Y = is given 
by 



Re^(A»l,y = 0) = L(i_JL 



(3.22) 



up to terms of order A 4 . For Y ^> 1 the approach to 
the limit A — > is Drude-like, 



Rea xx (X< 1,F> 1) 



e 2 2Fln2 
Il 2 + 4' 



(3.23) 



Second, we fix A = ui/rj and discuss the Y = l/{(3tvq) 
dependence of the real part of the longitudinal conduc- 
tivity from Eq. (|3 . 5|) . When Y ^> 1 as happens when the 
temperature /3 _1 is much larger than the energy resolu- 
tion hrj, {h/e 2 ) Re <r xx {X, Y) ~ gY for any fixed value of 
A where g is some constant. In particular, when A = 0, 
we find 

R C( i ;c;c (A = 0,y»l) = ^^y (3.24) 

to leading order in Y ^> L 23 i 44 In the opposite limit of 
Y <C 1, the conductivity approaches a finite value given 
by Eq. (pOTj) with Y = when A < 1, which is thus 
an increasing function of A <C 1 in agreement with Fig. 
QJb). The dependence on Y of Re a xx (A, Y) is monotonic 
increasing if A = while it is nonmonotonic for the finite 
values of A given in Fig. [Tfb). 
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FIG. 1: (Color online) Numerical integration of Eq. (|3.5[) 
when /i = v = x as a function of X = w/>) (a) and Y — 
{Ph-q)^ 1 (b). The dc limit u/rj — in (b) is obtained from 
numerical integration of Eq. (|3.17[) . 



IV. 



LANDAUER CONDUCTANCE 



In this section we are going to reproduce the calcula- 
tion of the Landauer conductance for a single massless 
Dirac fermion from Rcfs. iandi using the tools of CFT 
subjected to boundary conditions that preserve confor- 
mal invariance. Although the direct methods of Refs. [H 
and @ are both elegant and physically intuitive in the 
ballistic regime, we are hoping that the CFT approach 
might lend itself to a nonpcrturbative treatment of cer- 
tain types of disorder. 



Definition 



(4.1) 



where tr' denotes the trace over all channels in the 6th 
lead. The Landauer conductance (|4.1[) can be expressed 
in terms of the bilocal conductivity tensor £ of Eq. 
(JSH), according to2L4S.iL 



dS£^£^(r,r';e = 0,w = 0,»7 = 0), 

(4.2) 

where r (r') is constrained to lie on the interface between 
the ath (6th) lead and the sample, and J dS£ represents 
integration over the oriented interface between the sam- 
ple and the ath lead. 

To define the longitudinal Landauer conductance we 
choose for the sample the surface of a cylinder of length 
L x and of perimeter L y to which we attach two ideal 
leads l x and l 2 at the left end x L = —L x /2 and right end 
x R = +L x /2, respectively^ 

For the free Dirac Hamiltonian (|2.ip . the dimensionless 
conductance along the x direction, 



9 XX 



(h/e 2 )G L xx ee [h/e 2 )G\_ 



■2i 



(4.3a) 



can then be expressed in terms of the single-particle 
Green's function of Eq. 



9 XX 



where 



--(fovp) 2 <j> dy <j) dy' 

Ly o 

xtr [S* =0 (ry;0)a + S* =0 (r',r;0)<7_], 



(4.3b) 



(x L ,y), 



{x Rl y'), and er ± = a x ± ic 



We 

made use of the chiral symmetry and of ^^(r, r';e = 
0,lj = 0,n = 0) = S MM (r',r; £ = Q,w = 0,r) = 0). The 
single-particle Green's functions that enter Eq. (|4. 2|) arc 
obtained by solving the Schrodinger equation for the en- 
tire system, including the leads4^ For convenience, we 
assume that the leads also respect chira l 27 i 28 ' 29 i 30 sym- 
metry. The imaginary part Hrj of the energy can be set 
to zero in the sample, since the ideal leads broaden the 
energy levels in the sample. 

In the sequel, we will use Nazarov's technique to derive 
the following expressions for the dimensionless conduc- 
tance along the x direction, 



9 XX 



= ihL + 0[(L y /L x n 



(4.4) 



In order to define the Landauer conductance, we con- 
sider a finite region, the sample, described by the Dirac 
Hamiltonian, and attach a set of leads (or reservoirs) 
i x , l 2 , ■ ■ ■ to the sample. Propagation in the leads obeys 
different laws than in the sampled Then, the (dimcn- 
sionfull) conductance G^^ b for the transport from the 
ath to 6th lead is determined from the transmission ma- 



Thus, since the longitudinal conductivity a xx can be ex- 
tracted from the conductance in the anisotropic limit via 
9 xx = a xx L v /L x where L x L y , we recover from (|4.4[) 
the result (|3.20[) for the longitudinal Kubo dc conductiv- 
ity. 

The transverse Landauer conductance g xy can be de- 
fined by taking the sample to be a rectangular region 
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[-L x /2,+L x /2] x [-L y /2,+L y /2] and attaching four 
ideal leads to each edge. As is the case for the Kubo 
conductivity, we are going to show that 



(4.5) 



9 xv = °- 



We now turn to the derivations of Eqs. (|4.4p and (|4.5p 
for which we shall set 



h = v F = — e = 1 



(4.6) 



unless these constants are written explicitly. 



Nazarov formula 



Following Rcfs. IMij fl0lfli1fl2l . we introduce the generat- 
ing function for the transmission eigenvalue density 

7(B a Dct(l-7 £ 7^S^(0)^(0)) 

Here, Det refers to the functional determinant over all 
spatial coordinates (both inside and outside of the sam- 
ple) and spinor indices. We have also defined 



Vr/l = i^±S(x-x R/L ). 



(4.8) 



Finally, the source terms are parametrized by 9 F and B 
as 



1r = tan—, 



7l = sin cos 



(4.9) 



tan — — , Ct = sin — — cos — — . 
2 L 2 2 



The Landauer conductance is then given by 
dZ 



9 XX 



9{l L l R ) 



(4.10) 



Furthermore, if the transmission probability T n in chan- 
nel n (transmission eigenvalue), the nth positive real- 
valued eigenvalue of the product of transmission matrices 
entering Eq. (|4.ip in descending order, is written as 

T n =: cosh- 2 (0„/2), (4.11) 

then the density of transmission eigenvalues, 

n 

is given by 

p(9) =— [F(6 + i0 + + in) - F(9 - i0+ - in)] , 



F(0) 



_d_ 

del 



idO t 



Z{9 



F' U B, 



(4.12b) 



Once the density of transmission eigenvalues (|4. 12a|) is 
known, we can compute the Landauer conductance (see 
Ref.Hl) 



9xx ' 



the Fano factor 



F . 



En T n(l-TJ 



(4.13) 



(4.14) 



y t 

and other observables in terms of it. 

The essential step is to express the ratio of two determi- 
nants (|4.7p by fcrmionic and bosonic functional integrals 
as 

Z(6 F ,9 B ) = Z F (6 F )xZ B (6 B ), 
Z F = J® [x±,X±] e~ S ' F , Z B = Jv [ip ± ,tp ± ] 



where 



iS F = 



[X+,X- 



and 



(4.15a) 

ir){r) - 5i j R v R \ f X+ 

l L v L Mr) - H J V X- 

(4.15b) 



iri(r) ~ % CrVr _ 
( L v L ir](r) - 5i 



(4.15c) 

Here, f = J d 2 r denotes the space integral over the sam- 
ple and over the leads, {x a iXa) IS a P a i r °f two indepen- 
dent two-component fermionic fields, and {<p a ,tp a ) is a 
pair of two-component (complex) bosonic fields related 
by complex conjugation (cp a = tp* al a — ±) . In the func- 
tional integral, !K represents both the sample and leads, 
i.e., 3~C = IK inside the sample. Similarly, the smearing 
r](r) is zero in the sample but nonvanishing in the leads. 
We now turn to the modeling the leads. 



C. Boundary conditions 

There is quite some freedom in modeling the 'ideal' 
leads connected to the sample. In Ref. for example, 
propagation in the leads is governed by the nonrclativis- 
tic Schrodinger equation. In Ref. @ on the other hand, 
propagation in the leads is governed by the Dirac equa- 
tion with a large chemical potential. 

We are going to use this freedom to choose yet a third 
model for the leads. Wc demand that a Dirac fermion 
cannot exist as a coherently propagating mode in the 
leads. This can be achieved by choosing 

'K = 'K (4.16) 

in both the leads and the sample while using the smearing 
to distinguish between the sample and the leads, 



T)(r) 



0, r in the sample, 
oo, r in the leads. 



(4.17) 
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This choice for modelling the leads will be justified a 
posteriori once we recover from it the results of Refs. [H 
and0. (In the following we will use a cylindrical sample.) 

The spirit of the choice (|4.17p is similar to the prescrip- 
tion used in the nonlinear a model (NLcrM) description 
of weakly disordered conductors weakly coupled to ideal 
leads. 49 In the NLcrM for the matrix field Q. leads are 
represented by a boundary condition Q(x = ±L x /2) = A 
where A is a fixed matrix in the symmetric space of which 
Q is an element. The matrix Q describes the interact- 
ing diffusive modes of a weakly disordered metal. In a 
loose sense one may be able to think of this boundary 
condition as prohibiting coherent propagation of these 
diffusive modes in the leads. 

For a metallic sample (with a finite Fermi surface) in 
the ballistic regime that is weakly coupled to the leads, 
charge transport is strongly dependent on the nature of 
the contacts and the leads. On the other hand, for a 
metallic sample in the diffusive regime and not too large 
couplings to the leads, the conductance is mostly deter- 
mined by the disordered region itself. (See Ref. and 
references therein.) The conductivity of ballistic Dirac 
fermions in two dimensions is of order 1. Transport 
should thus behave in a way similar to that in a diffu- 
sive metals We would then expect that the microscopic 
modeling of the leads should have little effects on the 
conductance; i.e., the conductance should depend only 
on the intrinsic properties of the two-dimensional sample 
such as the conductivity. Reassuringly, it has been ob- 
served by Schomerus that transport in graphene is largely 
independent of the microscopic modeling of the leads 
Correspondingly, we will show that our model for the cou- 
pling between the sample and the reservoirs (|4T7l) leads 
to conformal invariant (i.e., scale-invariant and hence a 
renormalization group fixed point) boundary conditions 
to the supersymmetric field theory (|4.15|) . 



where 



The condition (|4.17[) suggests that the effects of the 
leads are equivalent to singling out a special configuration 
of the fields in the leads through the condition of a saddle 
point. To investigate the saddle-point condition implied 
by the leads (|4.17p . we introduce first the chiral base ip^, 
ip, j3', and (5 defined by 



{pi pi) = ^ x v^, ik) = 



(4.18) 



S 



(0) 



c'l 

Op 



CO 
F 



(4.19b) 



Q = ± 



\—iplip-6(x - x R ) + ^iplip + 5(x - x L ) 



and 



S 



(0) 



c"/ 
,5 B 



E / -{&w a +& a w a ), 
E f^Wa+PM, 



(4.19c) 



a=±"" 

^-I3\p_5{x - x R ) + ^Pl(3 + 5(x - x L ) 



The actions Sp^ and Sg^ give two copies of the Dirac 
fermion CFT (c = 1) and bosonic ghost CFT (c = — 1), 
respectively] 25 i 52 

In the leads, the field entering the functional integrals 
must then satisfy 

E + ^a + Pl'Pa + PlPa) = 0- (4.20) 

a=± 

Possible solutions to the saddle-point equations (|4.20[) are 



Pa = ±i0 a 



(4.21) 



Not all solutions (|4.21|) yield the desired Landauer con- 
ductance. One choice that does, as will be shown in Sees. 
IIVDI and IIVEI below, amounts to the boundary condi- 
tions 



tp a (x = TL x /2, y) = Ti4> a (x = TL x /2, y), 
P a (x = TL x /2,y) = T0 a (x = TL x /2,y), 



(4.22) 



with a = db. These boundary conditions break the fac- 
torization into a holomorphic and antiholomorphic sector 
present in the bulk. This is not to say that conformal in- 
variance is broken, however, as it is possible to eliminate 
one sector (say the antiholomorphic one) altogether in 
favor of the other (say holomorphic), thereby yielding a 
chiral conformal field theory.— 

At last, we need to impose antipcriodic boundary con- 
ditions in the (periodic) y direction of the cylinder, 



with a = ± in terms of which 



Z B = [ V[p ± ,P ± ]e- s - } - s -~ s -, 



(4.19a) 



fa 
fa 

Pa 
Pa 



fa 
fa 



(x,y + L ), 



(x,y) = -^ > j(x,y + L v ), 



(4.23) 



for a = ± and ~L x /2 < x < +L x /2 and < y < L . 
Our choice of antipcriodic boundary conditions for the 
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fcrmionic fields ip a and ip a is the natural one if the y di- 
rection is thought of as representing a "time" coordinate. 
The choice of periodic boundary conditions can be imple- 
mented at the price of introducing an additional operator 
in the conformal field theory. However, the Einstein con- 
ductivity does not depend on this choice of boundary 
conditions. 



D. 



Landauer conductance 



Before using the generating function (|4.7[) to compute 
directly the density of transmission eigenvalues (|4.12p . we 
compute the Landauer conductance (|4.10p as a warm-up. 
Insertion of Eq. ftlty into Eq. P~T0|) yields 



9xx = ^ <f <f dydy' (r) (r')) o 

L v 

(4.24) 

with r = (-L x /2,y) and r' = (+L x /2,y'). The ex- 
pectation value (•••)o is performed here with the ac- 
tion Sy + from Eqs. (|4.19|) supplemented with the 
boundary conditions (|4.22[) and (|4.23p . The four-fermion 
correlation function in Eq. (|4.24p can be expressed in 
terms of two-point correlation functions given by 



VJz, y)^ a , (0, 0)\ = ((3 a (x, y)0l (0, 0) 



with a, a' = ± and where^ 

(-1)' 



o (4.25a) 



5 (x,y,L x ,L y ) := 



m£t -f sinh f- (x + iy + 2mL x ) 

" (4.25b) 

After combining Eq. (|4.24j) with Eq. (|4.25p . one finds 



9 XX 



2^ cosh" 2 ( (2n + l) 



71 = 



.26) 



Each transmission eigenvalue 



T n :=cosh- 2 ( (2n+l)^ 
L y 



n = 0,1,2,- •• 



(4.27) 

is twofold degenerate. We shall see that this degeneracy 
originates from the two species (a = ±) when deducing 



(2n+ 1) 



27ri„ 



n = 0,1,2,- •• 



(4.28) 



directly from Eq. (|4TT2]) . 

The Landauer conductance (|4.26|) is a monotonic de- 
creasing function of L x /L y . When the sample is the 
surface of a long and narrow cylinder, L x /L y 1. the 
conductance is dominated by the contribution from the 



smallest transmission eigenvalue and decays exponen- 
tially fast with L x /L y ^> 1, 







-GnLJL. 



(4.29) 



In the opposite limit of a very short cylinder, L x / L y <C 1, 

oo 

l„. r d\ 
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27TL, 



cosh A 







[Ly/L x ) 



(4.30) 



7T L„ 



[Ly/L x ) 



We now turn to the computation of g xy . To this 
end, we take the sample to be a rectangular region 
[-L x /2, +L x /2] x [-L y /2, +L y /2] and attach ideal leads 
to each edge. Instead of the antiperiodic boundary condi- 
tion (|4.23[) , we must treat the boundary conditions along 
the y direction on equal footing with the boundary condi- 
tions along the x direction; i.e., we impose the boundary 
conditions 



^,0*. y = =F V 2 ) = Ttyafc. y = ^ V 2 )' 

a (x, y = t V 2 ) = Ti&fo V = TLJ2), 



(4.31) 



with a = ± together with the boundary condi- 
tions (|4.22[> . Equations (|3~6)) and (|4~Tj) , when applied 
to fl^, give 



I 

: 2n* 



+ L y/ 2 +LJ2 

J dy J dx 1 

-LJ2 -LJ2 



(4.32) 



where r = (-L x /2,y) and r' = (x',+L y /2). The ex- 
pectation value (•••)(} is performed here with the ac- 



tion Sp + Sg' supplemented with the boundary con- 
ditions (|4.22p and (|4.3ip . Using these boundary condi- 
tions, we can remove the right movers at the interfaces 
x = —L x /2 and y = +L y /2 with the result 



,(0) 



Qxy 



o. 



(4.33) 



E. Twisted partition functions 



We now go back to the direct calculation of the 
density of transmission eigenvalues, Eq. (|4. 12a|l . from 
Eq. (|4.12bp . We proceed in two steps. 

First, we perform a gauge transformation (defined 
in Appendix [B]) on the integration variables in the 
fcrmionic and bosonic path integrals, respectively, that 
diagonalizes the fermionic and bosonic actions 



S 



(0) 



:(o) 



C7 
° F 



S 

si 



F • 



(4.34) 
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In doing so the boundary conditions (|4.22[) that imple- 
ment the presence of the leads are changed to 

ip a (x = -L x /2,y) = -iTp a (x = -L x /2,y), 

^ a (x = +L x /2,y) = +i e + iae ^ a (x = +L x /2,y) 

(4.35a) 

and 

I3 a (x = -LJ2, y) = -ifrjx = -LJ2, y), 

a (x = +Lj2,y) = +ie +ae *P a (x = +L x /2,y), 

(4.35b) 

with a = ±, for the "gauge-transformed" fields. 

Second, we introduce the four independent partition 
functions Zp, Zp, Z B , and Z B describing two species 
(±) of fcrmionic (F) and bosonic (B) free fields (with 
holomorphic and antiholomorphic components) satisfy- 
ing the boundary conditions (|4.35| and (|4.23| . These are 
equivalent to four independent partition functions anti- 
^F:ch' Z F-,ch> Z B;dv an d Z B.ch describing free holomor- 
phic fields that fulfill the boundary conditions 



one verifies that 



ip ± (x,y + L y ) = - ip±(x,y), 



iP ± (x + 2L x ,y) 



^±{x,y), 



and 



(3±(x,y + L y ) = - /3 ± (x,y), 
P ± (x + 2L x ,y) = -e ±e °p ± (x,y), 



(4.36a) 



(4.36b) 



We have thus 



with — L x < x < L x and < y < L, 
traded the antiholomorphic sector in favor of a cylin- 
der twice as long and a change in the boundary condi- 
tions (|4.35[) implementing the presence of the leads. 

According to Ref. [HH, the chiral partition functions 
Zp. ch and Z|. ch are given by 



2nL r . 



p{6) = 2£*(0-(2n+l)-^). (4.39) 



The origin of the twofold degeneracy is the fact that the 
± species have decoupled. 

As it should be, the Landauer conductance is 
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d9p(9) cosh~ 2 (9/2 



2]Tcosh~ 2 ( (2n+l) 



L,, 



(4.40) 



Equation (|4.40|) implies that the density of transmis- 
sion eigenvalues is uniform, 



P(0) 



(4.41) 



in the limit nL x <C L y . In this sense, the transmis- 
sion eigenvalue density for the massless Dirac equation 
in a sample with the topology of a short cylinder agrees 
with that of a disordered metallic wire in the diffu- 
sive regime.— This is why transport for ballistic Dirac 
fermions is similar to mesoscopic transport in disordered 
quantum wires. 

The generating function technique can also be applied 
to g xy . If we follow the discussions for g xx , one simply 
finds that the partition function is actually independent 
of 6 FB . 



V. CHIRAL DISORDER WITH 
TIME-REVERSAL SYMMETRY 



zt ch = V 



n=0 



(4.37) 



n J ~ 



n=0 



+ e ±e Bq n +h l + e T e Bq n +V 



up to factors that cancel each other when we combine 
the fermionic and bosonic partition functions. We have 



introduced the variable q 



ZpZg = 1, separately for each species a = ±, when 
the boundary conditions are the same for fermions and 
bosons, i.e., when i9 F = 9 B , as it should be a conse- 
quence of global supcrsymmetry. With the help of 



^no+s 2 * 1 ) n 



n=0 



? I;ch=n( 1 + 

n=0 



n=0 




We devote this section to calculating the first-order 
correction to the Einstein conductivity induced by a weak 
real-valued random (white-noise) hopping amplitude be- 
tween nearest-neighbor sites of the honeycomb lattice at 
the band center. We are going to show that the Ein- 
stein conductivity is unchanged to this order. We then 
go on to show that the Einstein conductivity is an ana- 
lytic function of the disorder strength. 

We start from a single spinless fermion hopping be- 
tween nearest-neighbor sites of the honeycomb lattice at 
the band center. The hopping amplitudes are assumed 
real with small random fluctuations compared to their 
uniform mean. This model was introduced by Foster and 
Ludwig in Ref. [H. 

For weak disorder, this model can be simplified by lin- 
earizing the spectrum of the clean limit at the band cen- 
ter. In this approximation the clean spectrum is that 
of two flavors of Dirac fermions, each Dirac fermion en- 
coding the low-energy and long-wavelength description 
of the conduction band in the valley with a Fermi point. 
Weak disorder induces both intravalley and intervalley 
scatterings whose effects involve the two flavors of Dirac 
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fermions. The characteristic disorder strength for intra- 
node scattering is p A ; that for internode scattering is 
g M . 56 i 57 Without loss of generality, we shall set g A = 
and concentrate on g M > 0*££ In a fixed realization of 
the disorder, the single-particle Green's functions at the 
band center can be derived as correlation functions for 
fermionic (tp^,tp) and bosonic (ghost) (0^,0) variables 
from the partition function Z = Z F x Z B with 



Z F [m,m,A,A] = / D[ip^,ip] e"^"^, 



Z B [m,fh,A,A] = / D[0\0] 



-S R ^ R 



(5.1a) 



The low-energy effective (Dirac) action for the fermionic 
part is given by 

2 

S f = l^zZ + + rH?d + A)4, a 



"t , 



(5.1b) 



and 



S 



f = / ^ (V- 1 ^ 21 + - 1M1 - i> 2 Tpi) ■ (5-ic) 



The low-energy effective action for the bosonic part is 
given by the replacement S F —> S B under 

— > {p a \p a \(3 a j a ) (s.id) 

with a = 1, 2 and 



25. r_ Ait 

2tt 



The Abclian gauge fields A and A are source terms for the 
paramagnetic response function. The disorder is realized 
by the complex-valued random mass m and its complex 
conjugate m which obey the distribution law 

P\m\ oc / Dim, ml exp ( / mm) . (5. If) 

J y \ 2g M J r J 

The time-reversal symmetry of the lattice model has be- 
come 

-> -i/3 2 +, /3 2 t - +i/3 1 t, fl x -> +i(3 2 , /3 2 - -i/St, 

(5-lg) 

with the same transformation laws for the fields with 
overbars. The sublattice symmetry of the lattice model 
has become 

fl«t /3 »t j _ ^t, ^ ^ p a _ _^ 

(5.1h) 



for a = 1,2. 

We want to compute the first-order correction to the 
mean Einstein conductivity. To this end, we must double 
the number of integration variables in Eq. (|5.1ap . This 
is so because the response function is the product of two 
single-particle Green's functions. This is achieved by ex- 
tending the range a = 1, 2 of the flavor index in Eq. (|5 . 1 1) 
to a := (a a) with a = ± being a color index, one for 
each of the two single-particle Green's functions. [The 
same doubling of integration variables was introduced in 
Eq. (|4.15p .] We shall also use the more compact notation 
by which the capital latin index "A" replaces the original 
two- flavor indices a = 1, 2 in that it also carries a grade 
which is either when we want to refer to bosons - say, 
Y/a n = 0, i - or 1 when we want to refer to fermions - 
sa Y: X(Aa) = VVaa) ■ ^ * s tnc E T& dc of the indices A that 
enters expressions such as (— ) A . Correspondingly, we 
shall use the collective index 21 := (Aa) to treat bosons 
and fermions with the flavor index a = 1, 2 and the color 
index a = ± at once. 

Since we are only after the mean response function 
(and not higher moments), it can be obtained from 



Z[A,A] := J V[ X \x,x\x] exp(-S*-S v -S u ), 



S " : ~ I 2^°"' 



9m ,n 
2^ Um ' 



(5.2) 



Here, the interaction induced by integrating over the ran- 
dom mass with the probability distribution (|5.1fj) is de- 
scribed by 



Om :=X a X B+ X*X at (-l) A , 
while the smearing bilinear is 

Or, :=<*1,A<*1,B ( iT y) AB X^ + X» {^y) ^ X<B 

- S ,a s o,b X at {t x ) ab X St + Xa (^) AB X<b 



(5.3) 



(5.4) 



According to Eqs. (|3.17[) and (|3.18[) . the Einstein con- 
ductivity (|3.19p can be expressed in terms of the current- 
current correlation function — (ir 2 h 2 /e 2 )T, xx (r, 0). The 
latter function can be chosen to be represented in terms 
of bosonic variables, yielding 



J, 



'(1+) 



(1-) 

2-) 

'(2+) -r -(2+) 



(r) ( m + f ( n ) (o) 



(2+)\ 



J^ + J^Ur) (J'^ + Jfj (0) 



(5.5a) 
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[This is a generalization of Eq. (3.6) in which the single- 
particle Green's functions and currents are 4x4 matri- 
ces.] Here, wc have introduced the currents 

t'(oV) — a fl(o'a')t T'O'a') ._ 3 Ma'a')] 

J (aa) - — P(oo) P ' J (aa) '~ P(aa) P 

(5.5b) 

with the flavor indices a, a' = 1,2 and the color indices 
a, a' = ±, while the expectation value refers to 

((•■•)) := Jv[x\x,x\x] exp(-5 -5 M )(...) 

(5.5c) 



with S :— S t + SL. In the clean limit, the Einstein con- 
ductivity 

Be^^JJ^ry) (5.6) 

is given by twice the value of Eq. (|3.20|) at zero tempera- 
ture. This is understood as follows. The bilocal conduc- 
tivity reduces to computing the free-field expectation val- 
ues with the action S of bilincars in the normal-ordered 
current (|5.5b[) . By Wick's theorem, the bilocal conduc- 
tivity can be reduced to the products of pairs of free-field 
propagators. The relevant free-field propagators are 



(/ 3(2o)t (r) ^(W)t (0)) = (0(a«)t (r)/3 (ia')t (o)) = (f3 {la) (r)P {2al) (0)) = <% a) (r)/3 (2Q ,)(0)> = 5 aa ,2ir 



■k l V 



rj 2 + k 2 

(5.7a) 



and 



(P {aa) Hr)P {a , al) (0))=S aal S aa , 2m 



e^ k V-i, (P {aa} Hr)P (a , al) (0)) = S aa , S aa , 2m 



k V 2 + k 2 ' 



- ir - k x y , (5.7b) 



k V 2 + k 2 



with a, a' = 1,2 and a, a' = ± for the Einstein conduc- 
tivity. 

We turn next to the first-order correction in powers of 
g M of the mean Einstein conductivity and show that it 
vanishes. This is understood as follows. By translation 
invariance, the integration over r' in Eq. (|5 . 6[) yields 



Reosa. = ^ / S„(r,0), 



(5.8a) 



where the bilocal conductivity can be decomposed into 
three contributions, 

E xx (r, 0) = s£ h >(r, 0) + S^(r, 0) + ££ a >(r, 0). (5.8b) 

The holomorphic contribution is 



S£ h) (r,0) 



ir 2 H 2 ' 



BC'.JV^W^tO) ( 5 - 9a ) 



[summation convention over repeated indices a = (aa) is 
assumed on the right-hand side] with 



f> (!-)■ (2-) _ r> (2+), (1 + ) _ , 
n (l+), (2+) " ■°(2-), - 1 



(5.9b) 



the only nonvanishing coefficients. The antiholomorphic 
contribution is 



5t* a) M) 



TT 2 h 2 ~ 



B£'£(?£(r)f£(0)), (5.10a) 



with 



R (2~) _ o (2+), (1+) _ i 

(2+) - ^(2-), - 1 



(5.10b) 



the only nonvanishing coefficients. The mixed holomor- 
phic and antiholomorphic contribution is 

S £ a) M) = ^C£>*' (f£(r) f£(P)) , (5.Ha) 



with 



1 _ r (2-) _ „ (2-), 

1 (2+) ~ °(2+), (1+) 

_ r (2+). (1 + ) _ r (1 + ), (2+) 

_ °(2-), (1-) - °(1-), (2-) 



(5.11b) 



the only nonvanishing coefficients. The two-point func- 
tions (|5.9[) and (|5.10[) transform irreducibly and nontriv- 
ially under a rotation of the Euclidean plane r £ M 2 . 
Consequently, their separate contributions toEq. (|5.8p 
are vanishing and 

Re^ a = A^4fca)( r , ). (5.12) 

The first-order correction to the Einstein conductivity 
in the clean limit is 



<5Re< 



.9m e 



2ir 2 2ir 2 h 



^ a 2 , a 



2 r,Oh °1' a 3 



4^r)J: a /(0)O M (r') 



(5.13) 



Summation convention over repeated indices is assumed 
on the right-hand side. Carrying the double integration 
in Eq. (|5 . 13[) yields, with the help of Wick's theorem, 

SR ^=t^- h C ^^ 2 - (5 ' 14) 
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Summation convention over repeated indices is assumed 
on the right-hand side. Since C^'a" 1 = for any pair 
(a-^a-i) and (02^2)7 it follows that 



S Re a xx = 0. 



(5.15) 



Observe here that the factor tt 2 in Eq. (|5.14|) comes 
from the spatial integrations. It follows that the first- 
order correction to the Einstein conductivity is free from 
a logarithmic dependence on the ultraviolet cutoff. The 
correction of order g 2 ^ is also free from a logarithmic di- 
vergence but nonvanishingi^ These results are special 
cases of the fact that the Einstein conductivity must be 
an analytic function of the coupling constant g M . In- 
deed, it was shown in Ref. that the action + S M 
of Eq. (|5.2j) has the symmetry group GL(4|4) while the 
sector of the theory that carries no U(l) Abelian charge, 
the so-called PSL(4|4) sector, is a critical theory. One 
consequence of this is that the beta function for g M van- 
ishes to all orders in g M as M belongs to the PSL(4|4) 
sector. Another consequence is that the Einstein conduc- 
tivity must be an analytic function of g M as the bilocal 
conductivity also belongs to the PSL(4|4) sector.— 



VI. 



CONCLUSIONS 



We have shown how to compute the transmission 
eigenvalues for a single massless Dirac fermion propagat- 
ing freely in two dimensions within a two-dimensional 
conformal field theory description in the presence of 
twisted boundary conditions. We hope that this deriva- 
tion,which is complementary to the ones from Refs. [H 
and ly using direct methods of quantum mechanics, can 
be generalized to the presence of certain types of disorder 
so as to obtain nonpcrturbativc results. 



We have also shown that the Einstein conductivity, 
which is obtained from the regularization of the dc Kubo 
conductivity in terms of the four possible products of ad- 
vanced and retarded Green's functions by taking the dc 
limit before removing the smearing in the single-particle 
Green's functions, agrees with the conductivity deter- 
mined from the Landaucr formula. 

Finally, we noted that, as a consequence of the fixed 
point theory discussed in Ref. [H, the Einstein conduc- 
tivity is an analytic function of the strength g M of the 
disorder which preserves the sublattice symmetry of the 
random hopping model on the honeycomb lattice. More- 
over, the first-order correction in g M to the Einstein con- 
ductivity was shown to vanish. 
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APPENDIX A: NUMERICAL INTEGRATION 
OF EQ. (^5|) WHEN fi = v = x 



J 



To evaluate numerically Eq. (|3.5[) when fj, = v = x and for any finite 77 > it is useful to perform the integration 
over momenta in 



+ OO 



Rea^r^H— J ds [f f3 (e + /_) - /,(,)] J dkk ^_______-_ - ______ 



e — ihr] 



e + ihr] 



e + hw — ihr] 



e + huj + ihr] 



(s + hoj — ihr]) 2 — (hvFk) 2 (e + huj + ihi]) 2 — (hvFk)* 



This gives 
h 



1 1 



e2 ^a xx (X,Y)- 2nX 



da [f fi=1 (M/Y)-f fi=1 (p/Y)] 



X 



4a 



Iln__ 

Aa(a 2 + 1) ' X(X 2 + 4) j 2 n 1 + q 2 



X 



1 



4(a 2 + l) X 



— ) (arctanp— arctanq) 



2a 



1 



A 2 +4 2a 



— I (arctanp + arctang) 



where 



X 



1 

phi}' 



X 



hq ' 



p := a 



q := a + 



X 



(Al) 



(A2a) 



(A2b) 
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APPENDIX B: GAUGE TRANSFORMATION 

The "gauge field" 7 H / L can be removed by a suitable 

gauge transformation, tp — * U F (x)tp, ip^ — ► tp^U F (x)^ 1 , 
where U F (x) — V F W F (x). The gauge transformation 
W F {x) is position dependent, off-diagonal in the ± sec- 
tor, and given by 




lR°d${. x ~ x r) 



dx'A F (x') 



(Bl) 



where T x represents x ordering, and 



°d = 




2 



R/L 



2 




(B2) 



R/L 



Here, 



R/L 



denotes the holomorphic / antiholomorphic 



sector whereas ■ ■ - ± represents the ± sector. On the other 
hand, V F is position independent, diagonal in the holo- 
morphic and antiholomorphic sectors, and given by 



V F = 



a F 
a P 



R/L 



sin 6 F 



1 e 



1 



sin u F 



(B3) 



Since the action in the bulk is diagonal in the ± sector, V 
does not affect the action in the bulk, while it is chosen 
to diagonalize the boundary conditions at both ends of 
the cylinder. 

The "gauge field" Cr/l can be removed with the help 
of the gauge transformation, (3 — > V B W B (x)f3, (3' — > 
^Wgix)-^ 1 where 



A B (x) = 



C,l°J{ x - x l) 
W B (x) =T x exp I + J dx'A B {x') 



(R<7d S ( x ~ X R) 




(B4) 



and 



_ a B u 
Vb ~ I a f 



isin^s ismOs 

(B5) 
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